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The characterization of quantum polarization of light requires knowledge of all the moments of the Stokes 
variables, which are appropriately encoded in the multipole expansion of the density matrix. We look into the 
cumulative distribution of those multipoles and work out the corresponding extremal pure states. We find that 
SU(2) coherent states are maximal to any order whereas the converse case of minimal states (which can be seen 
as the most quantum ones) is investigated for a diverse range of the number of photons. Taking advantage of 
the Majorana representation, we recast the problem as that of distributing a number of points uniformly over the 
surface of the Poincare sphere. 

PACS numbers: 42.25.Ja, 42.50.Dv, 42.50.Ar, 42.50.Lc 


Introduction .— Stokes variables constitute an invaluable 
tool for assessing light polarization, both in the classical and 
quantum domains. They can be efficiently measured and lead 
to an elegant geometric representation, the Poincare sphere, 
which not only provides remarkable insights, but also greatly 
simplifies otherwise complex problems. 

Classical polarization is chiefly built on first-order moments 
of the Stokes parameters: states are pictured as points on 
the Poincare sphere (i.e., neglecting fluctuations altogether). 
Nowadays, however, there is a general agreement that a thor¬ 
ough understanding of the effects arising in the realm of the 
quantum world calls for an analysis of higher-order polariza¬ 
tion fluctuations Eta. In fact, this is what comes up in co¬ 
herence theory, where, in general, one needs a hierarchy of 
correlation functions to specify a field. 

Recently, we have laid the foundations for a systematic so¬ 
lution to this fundamental and longstanding question [ 8 -JO |. 
The backbone of our proposal is a multipole expansion of the 
density matrix, which naturally sorts successive moments of 
the Stokes variables. The dipole term, being just the first- 
order moment, renders the classical picture, while the other 
multipoles account for the fluctuations we wish to scrutinize. 
Consequently, the cumulative distribution for these multipoles 
yields complete information about the polarization properties. 

This Communication represents a substantial step ahead in 
this program, as we elaborate on the extremal states for the 
aforementioned multipole distribution. We find that the SU(2) 
coherent states maximize it to any order, so they are the most 
polarized allowed by quantum theory. We determine as well 
the states that kill the cumulative distribution up to a given or¬ 
der M : they serve precisely as the opposite of SU(2) coherent 
states and hence can be considered as the kings of quantum¬ 
ness. Furthermore, employing the striking advantages of the 
Majorana representation m , these kings appear naturally re¬ 
lated to the problem of distributing A points on the Poincare 
sphere in the “most symmetric” fashion, a problem with a long 
history and many different solutions depending on the cost 
function one tries to optimize EHH. 


Polarization multipoles .— The quantum Stokes operators 
are defined as m 

S x = ^(a\d- +afa +), S y = ^(a + a^_ — d+d -), 

(1) 

S z = \(d+d+ — afd -), 

together with the total photon number A = a\a+ f-aLd-. 
Here, a+ and d- represent the amplitudes in two circularly- 
polarized orthogonal modes. We have that [a^d\] = <5^, 
k,£ G {+,—}, with h = 1 throughout and the superscript t 
stands for the Hermitian conjugate. The definition ([T]) dif¬ 
fers by a factor 1/2 from its classical counterpart [15], but 
in this way the components of the vector S = (S x ,S y ,S z ) sat¬ 
isfy the su(2) commutation relations: [5 x ,Sy] = iS z and cyclic 
permutations. For an A-photon state, S 2 = S(S-\- 1)11, where 
S = A/2, so the number of photons fixes the effective spin. 

Put in a different way, Q is nothing but the Schwinger rep¬ 
resentation of the su(2) algebra. Consequently, the ideas to be 
explored here are by no means restricted to polarization, but 
concern numerous instances wherein su(2) is the fundamental 
symmetry he 

In our case, [A,S] = 0, so each subspace with a fixed 
number of photons ought to be addressed separately. To 
bring out this fact more prominently, instead of the Fock 
states {|/z + , zz_)}, we employ the relabeling \S,m) = \n+ = 
S + m,n- = S — m), which can be thought of as the common 
eigenstates of S 2 and S z . For each fixed S , m runs from —S to 
S , and the states {| S,m)} span a (2 S + 1)-dimensional invari¬ 
ant subspace m. 

As a result, the only accessible information from any den¬ 
sity matrix p is its block-diagonal form p po i = ©y P^, where 
p ® is the density matrix in the subspace of spin S. This p po i 
has been termed the polarization sector m or the polariza¬ 
tion density matrix he It is advantageous to expand each 
p® as 

2s k 

p (s, = E E 

K=()q=-K 
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rather than using directly the basis {| 5,ra)}. The irreducible 
tensor operators tjfj are I201I2F 1 
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with Cg^ K being Clebsch-Gordan coefficients (0 < K < 25). 
These tensors form an orthonormal basis and have the right 
properties under SU(2) transformations. The crucial point is 

yv 

that w can be jotted down in terms of the ^th power of the 
Stokes operators. 

The expansion coefficients = Tr[p^ are known 

as state multipoles. The quantity Y, q \ pKq 1 2 £ au S es the state 
overlapping with the Kth multipole pattern. For most states, 
only a limited number of multipoles play a substantive role 
and the rest of them have an exceedingly small contribution. 
Therefore, it seems more convenient to look at the cumulative 
distribution 

<’=I E lpgl 2 > (4) 

K=1 q=—K 


which sums polarization information up to order M (1 < M < 
25). Note that the monopole K = 0 is omitted, as it is just a 
constant term. As with any cumulative distribution, g/^ ’ is a 

monotonically nondecreasing function of the multipole order. 

(s) 

Maximal states .— The distribution g/^ ’ can be regarded 

( 5 ) 

as a nonlinear functional of the density matrix p . On that 

(s) 

account, one can try to ascertain the states that maximize g/^ ’ 
for each order M. We shall be considering only pure states, 
which we expand as |T*) = |5,m), with coefficients 

T f m = (5,m|T / ). We easily get 
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The details of the calculation are presented in the Supplemen¬ 
tal Material. We content ourselves with the final result: the 
maximum value is 

( 5 ) 25 [r(25+l)] 2 

M 25+1 T(2S-M)T(2S + M + 2) ’ V 

and this happens for the state |5,±5), irrespective of M. 
Since g/^ J is invariant under polarization transformations, all 
the displaced versions |0,0) = (1 + |a| 2 ) _lS exp(a5 + )|5, —5) 
[with S± = S x ± iS y and the stereographic projection a = 

tan(0/2)e - ^] also maximize gA^\ In other words, SU(2) 

coherent states 10,0) l22l maximize gfjp for all orders M. 

It will be useful in the following to exploit the Majorana 
representation CD, which maps every (25+ 1)-dimensional 
pure state |*F) into the polynomial 


¥(«)= £ , 

m=—S 
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(5 —m)!(5 + m)! 
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Up to a global unphysical factor, I^P) is determined by the set 
{(Xi} of the 25 complex zeros of 'P(oc), suitably completed 
by points at infinity if the degree of 'P(a) is less than 25. A 
nice geometrical representation of IT 1 ) by 25 points on the 
unit sphere (often called the constellation) is obtained by an 
inverse stereographic map of {«/} i-+ {0/,0/}. For SU(2) co¬ 
herent states, the Majorana constellation collapses to a sin¬ 
gle point. States with the same Majorana constellation, irre¬ 
spective of its relative orientation, share the same polarization 
properties. 

TheSU(2) Q-function, defined as Q(0, 0) = | (0,0 (T 1 ) | 2 , is 
an alternative way to depict the state. Although 2(0,0) can be 
expressed in terms of the Majorana polynomial [and so {(+■} 
are also the zeros of Q(0, 0)], sometimes the symmetry group 
of IT 1 ) can be better appreciated with this function, which can 
be very valuable. 

Minimal states. — Next, we concentrate on minimiz¬ 
ing gfjjp. Obviously, the maximally mixed state p® = 
2 ^yll 2 S+i kills all the multipoles and so indeed causes to 
vanish for all M, being fully unpolarized l23l[24j . Nonethe¬ 
less, we are interested in pure Mth-order unpolarized states. 
The strategy we adopt is thus very simple to state: starting 
from a set of unknown normalized state amplitudes in Eq. ([5]), 

which we write as T f m = a m + ib m (a m , b m E M), we try to get 

(s) 

g/^ — 0 for the highest possible M. This yields a system 
of polynomial equations of degree two for a m and b m , which 
we solve using Grobner bases implemented in the computer 
algebra system MAGMA ll25ll . In this way, we get exact alge¬ 
braic expressions and we can detect when there is no feasible 
solution. 

Table [T] lists the resulting states (which, in some cases, are 
not unique) for different selected values of 5 l26l . We also in¬ 
dicate the associated Majorana constellations. For complete¬ 
ness, in Fig. [T] we also plot the constellations as well as the 
2-functions for some of these states. 

Intuitively, one would expect that these constellations 
should have the points as symmetrically placed on the unit 
sphere as possible. This fits well with the notion of states of 
maximal Wehrl-Lieb entropy ED. In more precise mathe¬ 
matical terms, such points may be generated via optimization 
with respect to a suitable criterion d. Here, we explore the 
connection with spherical ^-designs l28l . which are patterns 
of N points on a sphere such that every polynomial of degree 
at most t has the same average over the N points as over the 
sphere. Thus, the N points mimic a flat distribution up to order 
t , which obviously implies a fairly symmetric distribution. 

For a given 5, the maximal order of M for which we can 

(s) 

cancel out g’ does not follow a clear pattern. The numer¬ 
ical evidence suggests that M max coincides with t m ax in the 
corresponding spherical design, but further work is needed to 
support this conjecture. 

The simplest nontrivial example are two-photon states, 5 = 
1. We find only first-order unpolarized states: these are bipho¬ 
tons generated in spontaneous parametric down-conversion, 
which were the first known to have hidden polarization f30l . 
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TABLE I. States that kill ; for the indicated values of 5. In the second column, we indicate the order M, which we conjecture is the highest 
possible. We give the nonzero state components l P m (m = —S ,... ,5) and the Majorana constellation. We include the associated spherical t- 
design (with the maximal t value) and the Queens of quantumness (with their unpolarization degree), “same” “similar” and “different” always 
refer to the closest description column to the left. 


5 

M 

State 

Constellation 

Design 

t 

Queens 

M 

1 

1 

¥0 = 1 

radial line 

same 

1 

same 

1 

3 

2 

1 

*-i=*i = A 

equatorial triangle 

same 

1 

same 

1 

2 

2 

(Nlm 

II 

<N 

-1^ 

II 

7 

tetrahedron 

same 

2 

same 

2 

5 

2 

1 

*-i=*i = 7S 

equatorial triangle + poles 

same 

1 

same 

1 

3 

3 

^-2-^ = 4 

octahedron 

same 

3 

same 

3 

7 

2 

2 

1 

II 

II 

II 

l^toj 

two triangles + pole 

similar 

2 

equatorial pentagon + poles 

1 

4 

3 

^-4 = ^4 = VX *0 = V7 

cube 

same 

3 

see lf29l 

1 

9 

2 

2 

v 9 =¥<> = 4 , m 3 =^3 = 4 

2 2 V 6 2 2 V5 

three triangles 

similar 

2 

similar 

1 

5 

3 

= 75, *0 = 73 

pentagonal prism 

similar 

3 two staggered squares + poles 

1 

11 

2 

3 

1 

HE 

II 

HE 

II 

l 

II 

II 

pentagon + two triangles 

similar 

3 

similar 

1 

6 

5 

•P-S = -^5 = #, 'i '0 = -4 I 

icosahedron 

same 

5 

same 

5 


7 4 *P _ 6 = ¥ 6 = Up T / _3 = ¥3 = + iyj three squares + poles different 4 


xu _ / 12561757 ; / 512603 

Y 0 ~ V 163053000 W 2013000 

10 5 *F io = ^10 = a/ 7g73, *¥-5 = —^5 = a/§§§, 'Pq = \T^§h deformed dodecahedron similar 5 


With three photons, 5 = 3/2, we have again only first-order 
unpolarized states: the constellation is an equilateral triangle 
inscribed in a great circle, which can be taken as the equator. 
This coincides with the three-point spherical 1-design. 

For 5 = 2, the Majorana constellation is a regular tetrahe¬ 
dron: it is the least-excited second-order unpolarized state. It 
is not surprising that the tetrahedron is the 2-design with the 
lowest number of points. 

The case 5 = 5/2 does not admit a high degree of spherical 
symmetry: only first-order unpolarized states exist. There are 
neither five-photon M = 2 unpolarized states l3Tll32ll nor five- 
point 2-designs ED- 

When increasing the number of photons to six, 5 = 3, an¬ 
other Platonic solid appears: the regular octahedron. Now, we 
have the least-excited third-order unpolarized states, which, in 
addition, take on the maximum sum of the Stokes variances. 

For 5 = 7/2, an M = 2 constellation consists of the north 
pole, an equilateral triangle inscribed at the z = 0.2424 plane, 
and another equilateral triangle, with the same orientation 
(e.g., one vertex on the x-axis) at the z = —0.5816 plane. The 
spherical Ldesign has a larger separation between the trian¬ 
gles; but the corresponding Stokes vector does not vanish, so 
the Ldesign does not coincide with any unpolarized state. 

The next Platonic solid, the cube, appears when 5 = 4. The 
state is third-order unpolarized and its Majorana constellation 
coincides with the eight-point spherical 3-design, which is the 
tightest for this number of points. 

A nine-photon second-order unpolarized state, 5 = 9/2, is 
generated by three equilateral triangles with the same orienta¬ 


tion inscribed in the equator and in two symmetric rings. The 
highest nine-point spherical Ldesign has t = 2 and a similar, 
but not identical, configuration because the two smaller trian¬ 
gles are displaced by a larger distance from the equator than 
the previous constellation. As a consequence, the nine-point 
spherical 2-design is only first-order unpolarized. 

The Majorana constellation for a maximally unpolarized 
10-photon state (5 = 5) is similar to the matching spherical 
Ldesign and consists of two identical regular pentagons in¬ 
scribed in rings symmetrically displaced from the equator. 
The maximally unpolarized state has the two pentagons a 
bit closer to the equator than the spherical 3-design (that has 
M = 1). 

For larger photon numbers, the computational complexity 
of finding optimal designs becomes a real hurdle. The t- 
designs have been investigated in the range 2-100 and nu¬ 
merical evidence suggests that the optimal designs (in some 
instances, they are not unique) have been found f34l . How¬ 
ever, for some dimensions, e.g., 12 (5 = 6) and 20 (5 = 10), 
one would naively guess that the optimal designs fit with the 
icosahedron and the dodecahedron. For 5 = 6 this turns out 
to be a correct guess, the corresponding state is unpolarized to 
the same order as the spherical 5-design formed by the icosa¬ 
hedron. For 5 = 10 this intuition fails: the optimal Ldesign 
is indeed a dodecahedron, but this Majorana constellation is 
third-order unpolarized, whereas this is a spherical 5-design. 
If the dodecahedron is stretched (i.e., the four pentagonal rings 
that define its vertices are displaced against the pole), one can 
find a 20-photon fifth-order unpolarized state. 
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FIG. 1. (Color online) Density plots of the SU(2) Q functions for the optimal states in Table I for the cases S = 5/2,3,7/2,9/2,5, and 7 (from 
left to right, blue indicates the zero values and red maximal ones). On top, we sketch the Majorana constellation for each of them. 


To check the correspondence between unpolarized states 
and optimal ^-designs we look at dimension 14, which is the 
smallest number of points for which a spherical 4-design, but 
not a 5-design, exists. This consists of four equilateral tri¬ 
angles that are pairwise similar in size, displaced from the 
equator by the same distance, and rotated an angle ±a or ±/3 
around their surface normal, plus the two poles. The ^-design 
state is only first-order unpolarized, but if the spacing and tri¬ 
angle orientation is optimized, the design can be made third- 
order unpolarized. There is indeed a 14-photon state that is 
fourth-order unpolarized: its Majorana constellation is made 
of three quadrangles and the poles, but this is only a 1-design. 

To round up, it is worth commenting on the connections 
that our theory shares with two recently introduced notions: 
anticoherent states (351 and queens of quantumness (29) . For 
completeness, in Table [I] we have also listed the configura¬ 
tions and the degree of unpolarization for these queens. Anti¬ 
coherent states are in a sense “the opposite” of SU(2) coher¬ 
ent states: while the latter correspond as nearly as possible to 
a classical spin vector pointing in a given direction, the for¬ 
mer “point nowhere”, i.e., the average Stokes vector vanishes 
and the fluctuations up to order M are isotropic. The queens 
of quantumness are the most distant states (in the sense of a 
Hilbert-Schmidt distance) to the classical ones (states than can 
be written as a convex sum of projectors onto coherent states). 
In particular low-dimensional cases, these two instances co¬ 
incide with our optimal states. However, we stress that our 
theory is built from first principles, starting from magnitudes 
that are routinely determined in the lab. Besides, we have an 
algebraic criterion, namely, the vanishing of the cumulative 
multipole distribution, that can be handled in a clear and com¬ 
pact manner. 

When we interpret our (2S + 1)-dimensional subspace as 
the symmetric subspace of a system of S qubits, the kings ap¬ 
pear also closely linked to other intriguing problems, such 
as maximally entangled symmetric states (361 ED and k- 
maximally mixed states [38] 3.91. 

Applications .— The main goal of quantum metrology is to 


measure a physical magnitude with surprising precision by ex¬ 
ploiting quantum resources. In particular, tailoring polariza¬ 
tion states to better detect SU(2) rotations is quite a relevant 
problem with direct applications to magnetometry, polarime- 
try, and metrology, in general m. 

In this respect, MXW states [defined as \N00N) = (| S,S) — 

| S,—S))/y/2\ are known to be maximally sensitive to small 
phase shifts (i.e., to small rotations about the S z - axis) for a 
fixed excitation S m . This can be easily understood by look¬ 
ing at their Majorana constellation, which consists in just 2 S 
equidistantly placed points around the Poincare sphere equa¬ 
tor. Since a rotation around the S z axis is described by the uni¬ 
tary operator U(&) = exp(— i&S z /2), the states |N00N) and 
U(&) |NOON) are orthogonal for n/(2S). However, to make 
optimal use of a NOON state it is essential to know the rota¬ 
tion axis so as to ensure that the state is aligned with the axis 
to achieve its best sensitivity: the rotation resolution is thus 
highly directional for a NOON state. 

This is precisely the advantage of maximally unpolarized 
states: having a high degree of spherical symmetry, they re¬ 
solve rotations around any axis approximately equally well. 
This has been confirmed for the Platonic solids EB: Pla- 
tonic states saturate the optimal average sensitivity to rotations 
about any axis; NOON states outperform these states about one 
specific axis (42). Indeed, for the Platonic solids, rotations 
around all the facets normal axes map the Majorana constella¬ 
tion onto itself for rotations of 2;r/3 (tetrahedron, octahedron 
and icosahedron), ;r/2 (cube), or 27 t/ 5 (dodecahedron). It is 
clear that for other constellations and other rotation axes the 
Majorana constellation will only become approximately iden¬ 
tical, but the statement is more likely to hold true. 

In a different vein, we draw attention to the structural sim¬ 
ilarity between the kings of quantumness and quantum error 
correcting codes: in both cases, low-order terms in the expan¬ 
sion of the density matrices are required to vanish. 

As a final but relevant remark, we stress that all the ba¬ 
sic tools needed for our treatment (Schwinger representation, 
multipole expansion and constellations) have been extended 
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in a direct way to other symmetries, such as SU(3) US or 
Heisenberg-Weyl ED. Therefore, the notion of kings of quan¬ 
tumness can be easily developed for other systems. Work 
along these lines is already in progress in our group. 

Concluding remarks. — In short, we have consistently 
reaped the benefits of the cumulative distribution of polariza¬ 
tion multipoles, which is a sensible and experimentally real¬ 
izable quantity. We have proven that SU(2) coherent states 
maximize that quantity to all orders: in this way, they mani¬ 
fest their classical virtues. Their opposite counterparts, mini¬ 
mizing that quantity, are certainly the kings of quantumness 

Apart from their indisputable geometrical beauty, there 
surely is plenty of room for the application of these states, 
whose generation has started to be seriously considered in sev¬ 
eral groups. 
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(Grant Q-ESSENCE), and the Program UCM-Banco San¬ 
tander (Grant GR3/14) is gratefully acknowledged. GB thanks 
the MPL for hosting him and the Wenner-Gren Foundation for 
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Appendix: Optimal states. — We have to maximize the cu¬ 
mulative multipole distribution © for a pure state |*P) — 
Tsj n =-S l^m), which takes the form <0- If we use integral 
representation for the product of two Clebsch-Gordan coeffi¬ 
cients 03 , we get 


^{S) _ y y 2S + 1 y y 2K +\ ;< j,* ^ ^* 

K ~ L L Sn 2 L L 2S +I* m 

m,m =—S n,n =—S K=\)q=—K 

x jdRDl n {R)D s :, n ,{R)D K qq {R), (8) 


where D s mn are the Wigner D-functions and R refers to the 
three Euler angles (a, j3, y) and the integration is on the group 
manifold 


/ pin p7l pin 

dRf{R) = da <7/3 sin/3 / dyf{a, J3,y) . (9) 

7o 7o 7o 

Since 

D« q (R) = XkM, (10) 

q——K 


where Xk{®) is a SU(2) generalized character and 
cos(to/2) = cos(/3/2)cos[(a + y)/2], we rewrite g/jp as 

= £ ^J dRXK(co) | (^\T g s \W) | 2 . ( 11 ) 


and f g is the group action. Then, we observe that the above is 


4 S) = Tr 




M 


K =0 


2K +1 
8tt 2 


J dRXK{co)Tg 0: 


( 12 ) 


with 


1*)= i (-im m \S,rn). (13) 

m=—S 


Because the integral 

A* J dR XK((0)T g s ® tp = c K n K , (14) 

where IT^ is the identity on the (2 K + 1)-dimensional irre¬ 
ducible SU(2) subspace which appear in the tensor product of 
[i.e., Tr(n^) = 2K +1], then 

M 

^m= £<^|<^|n JC |^)|'l'). (15) 

K= 1 

Such overlap is maximized (all coefficients are the same) 
whenever in every subspace of dim 2K + 1 there is only one 
element from the decomposition IT 1 )IT 1 ), which is consistent 
with ( p~3] ). The only states that at decomposition on represen¬ 
tations produce a single state in each invariant subspace are 
the basis states |S,m), so that |*P) = (— 1 ) m \S, —m), then 



k= 1 


2K +1 
25+1 


r s-s 

^SS,K0 


(16) 


Since the maximum value of Cf® s _ m is Cf$ s _ s , the states 

fc\ 

| S, ±S ) maximize as heralded before. 
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